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In order to determine types of the orbital ordering in man- 
ganites, we study theoretically the polarization dependence 
of the anomalous X-ray scattering which is caused by the 
anisotropy of the scattering factor. The general formulae of 
the scattering intensity in the experimental optical system is 
derived and the atomic scattering factor is calculated in the 
microscopic electronic model. By using the results, the X-ray 
scattering intensity in several types of the orbital ordering is 
numerically calculated as a function of azimuthal and ana- 
lyzer angles. 
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I. INTRODUCTION 

Perovskite manganites and their related compounds, 
(Ri_ x A x )i+ n Mn n 03n+i (n = 1, 2, oo, R = rare earth ion, 
A — alkaline earth ion), have been studied extensively 
from the fundamental scientific view point as well as the 
technological one, sinc&,the discovery of the colossal mag- 
netoresistance (CMR).lil □ In order to explain the unique 
magnetic and transport properties, special attention is 
paid to the orbital degrees of freedom in Mn ions.ETtll 
The electron configuration of a Mn 3+ ion is represented 
by (^2 g ) 3 (e g ) 1 due to the strong Hund coupling. Since 
an electron occupies the two degenerate e g orbitals, the 
ion has the orbital degrees of freedom as well as spin 
and charge. It was supposed that the orbital rffidpr.- 
ing is realized in Insulating compounds RMnDsSoEjO, 
Ro.5Ao. 5 MnOja<BHo and Lao.5Sr1.5MnO4.EEl Further- 
more, in the metallic phase, roles of the orbital degrees 
in the magneiii^-pptical and transport phenomena were 
also stressedU'EZrO However, the direct observation of 
the orbital degrees was restricted experimentally to the 
case where the polarized neutron diffraction is utilized.Ed 

Recently, Murakami et al. have applied the anoma- 
lous X-ray scattering in order to detect the orbital o&, 
dering in single layered manganites Lao.5Sr1.5MnO4.Ell 
They focused on a reflection at (| | 0) point and ob- 
served a resonant-like peak near the K-edge of a Mn 3+ 
ion below about 200K. They further observed the unique 
polarization dependence which is attributed to the ten- 
sor character of the anomalous scattering factor. Since 
all Mn 3+ ions are equivalent, the reflection at (| | 0) 
is forbidden. Therefore, an appearance of the intensity 
implies that two kinds of orbital are alternately aligned 
in the Mn02 plane. This type of the orbital ordering 



is termed the antiferro-type orbital ordering, hereafter. 
The experimental results also imply that the dipole tran- 
sition between Mn Is and Mn 4p orbitals causes the 
scattering. The experimental, method was extended to 
Lai_ x Sr x Mn0 3 with x = O.OB and 0.12.E3 

In this paper, we study theoretically the polarization 
dependence of the anomalous X-ray scattering as a probe 
to identify types of the orbital ordering in manganites. 
Although the phcnomcnological model for the anoma- 
lous scattering factor was used to analyze the experimen- 
tal results of (§ | 0) reflection, it was not determined 
which type of the orbital ordering, (3d 3x 2_ r 2 , 3c?3j,2_ r gA 
or (3d z 2_ x 2,3d y 2_ z 2), is realized in Lao.5Sr1.5MnO4.EiJ 
As for the theoretical side, several types of the orbital 
ordering is predict ed. in Xhc wide range of the carrier 
concentrationBLniiO'Ota Once a method to identify 
types of the orbital ordering is established, this becomes 
applicable to the study of the orbital structure in not 
only manganites but also a wide range of transition metal 
oxides. Considering the fact that types of the orbital or- 
dering directly reflects on the anisotropgi_pf tensor ele- 
ment of the anomalous scattering factor JEjI 3 !] we investi- 
gate the polarization dependence of the scattering inten- 
sity. Through the numerical calculation in several orbital 
ordered cases, we propose the method to identify the or- 
derings. 

In Sect. II, the general formulae of the polarization 
dependence of the scattering intensity is derived. In Sect. 
Ill, the anomalous part of the-,atomic structure factor is 
calculated in a MnOe cluster J 3 -!! In Sect IV, the numerical 
results of the polarization dependence in several orbital 
ordered cases, which are calculated by using the results 
obtained in Sects. II and III, are presented. Sect. V is 
devoted to the summary and discussion. 

II. GENERAL FORMULAE 

In this section, we derive the general formulae of the 
polarization dependence of the scattering intensity in or- 
bital ordered manganites. The conventional experimen- 
tal arrangement, where the poiajfiaation measurement is 
carried out, is shown in Fig.l.E3E3 It mainly consists of 
a sample crystal and the polarization analyzer including 
analyzer crystal and photon detector. The polarization 
scan is characterized by two rotation angles, i.e., the az- 
imuthal angle (tp) and the analyzer one (if a)- The former 
is the rotation angle of the sample around the scattering 
vector K = k" — k', and the latter is that of the analyzer 
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around an axis which is parallel to the scattered photon 
beam. The incident beam from the synchrotron source 
is almost perfectly polarized in the horizontal plane, i.e., 
cr-polarization. The direction of the electric vector of the 
incident photon with respect to the crystalline axis is 
changed by the azimuthal rotation. The scattered pho- 
ton has both n- and cr-polarization components, due to 
the tensor character of the anomalous scattering factor, 
which are separated by the analyzer scan. 

Jjjjt&is optical system, the scattering intensity is given 



by! 



I(k>,k>',<p,v A )=J2\J2 M x^ A ^ k '> k ">^r > (!) 

A' A 

for A (A') =<7 or tt. M\>\(pa) is the scattering matrix in 
the analyzer and described by the analyzer angle and the 
scattering one 9 a in the analyzer as follows, 



M x <\{<pa) = F a 



COS if A 

sin ipA cos 29a 



— sin ifA 
cos (fA cos 29a 



(2) 



where Fa is the scattering factor in the analyzer crystal. 
In this paper, the scattering angle is fixed at 9a = tt/4, 
as chosen in the conventional experiments. The a- and 
7r-polarized components in the scattered photon are de- 
tected by the photon detector with <pa = and tt/2, 
respectively. A\ny(k", k 1 , ip) in Eq. ([!]) is the scattering 
amplitude with the incident (scattered) photon momen- 
tum k' {k") and the polarization A' (A") as a fuction of 
the azimuthal angle. It is given by, 



THY 1 ^ 



A"a" 



U{(p)VF(k',k>')V' i U{(py 



(COS ip 
simp 





(6) 



and the unitary matrix V describes the transformation 
from the coordinate in the crystallographic axis (a, b, c) to 
that in the laboratory system (ei,e*2,e3). The structure 
factor F a p(k' , k") in Eq.(^) is represented by a sum of 
the normal and anomalous parts of the scattering factor 
as follows: 

F aP (k>, k") =N^ e. i(g - kl ' yn fi«p(k>, k") , (7) 



i^cell 



-A'a' > 



and 



-Aa 



-Aa 



1 

sinf 





cos( 



(4) 



(5) 



with 



MW, k") = foi(k', k")S a0 + Af iaf3 {k>, k") . (8) 

Here, foi(k',k") and Afi a p(k',k") are the normal and 
anomalous parts of the atomic scattering factor, respec- 
tively, of the i-th atom defined in the (a, b, c)-coordinate. 
N is the number of the unit cell. When we determine 
the unitary matrix V and the atomic scattering factors, 
the scattering intensity is calculated as a function of the 
azimuthal angle ip and the analyzer angle ipA- 

It is instructive to demonstrate the polarization de- 
pendence of the scattering intensity in the following 
two cases: (i) (e"i,e2,e3) = (a,b,c) (azimuthal rota- 

' (a- 
h- 

axis). When the crystallographic axes coincide with the 
principle ones in the MnOg octahedron, the anomalous 
part of the scattering factor in Eq. (||) is diagonal as 
Afi a p{k',k") = 5 a pAfi aa (k',k"). In this condition, the 
scattering intensity is given by 



tion around the c-axis) and (ii) (ex,e2,e3) = (c, — ^ 
6), -^(a + b)) (azimuthal rotation around the a - 



(3) 



I(k>,k»,<p,ip A ) 



mc 
x \D„„ 



COS <p>A 



D 



tt(t sin p>A\ 



where F a p(k', k") is the structure factor of the sample 
crystal defined in the coordinate of the crystallographic 
axis described as (a,b,c). e^ a and ejj^ are the polariza- 
tion vectors of the incident and scattered photons, re- 
spectively, and their explicit forms are given by 



with 



and 



D aa = F xx cos ip + F yy sin ip 



(9) 



(10) 



= ( F xx + Fyy) sin ip cos ip sin 9 , (11) 
for case (i), and 



Daa = F zz cos 2 tp + -(F xx + F yy ) sin 2 <p , (12) 



and 



where the suffix a denotes the Cartesian coordinate in 
the laboratory system, i.e., (ei,e2,e3) shown in Fig. 1. 

e \a ^ s ^ ne transposed matrix of . The unitary matrix 
U(p) in Eq. (||) describes the azimuthal rotation of the 
sample around the e3-axis defined by 



Dna = 2 (Fxx + F ljy - 2F ZZ ) sin tp cos tp sin I 
+ o {—Fxx + Fyy) sin p cos 9 , 



(13) 



for case (ii). 
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III. ATOMIC SCATTERING FACTOR 

The normal and anomalous parts of the atomic scat- 
tering factor in Eq.(|^) are obtained by the perturbation 
calculation with respect to the electron-photon interac- 
tion. The normal part is given by the Fourier trans- 
form of the charge density pi in the i-th atom defined by 
foi = <J\fH{K = k" - #)|0) , where |0) (|/)) is the initial 
(final) electronic state with energy e The anoma- 

lous one is derived by the interaction between electronic 
current and photon and is expressed as follows oa 



Af ia p(k',k") 



m^l m a (~k>)\l)(l\j i0 (k")\O) 



e - ei - u k n - i5 

(f\Mk^)\l)(l\jU-k')\0) \ ^ 
e - £; + LOk' ~i5 f 



(14) 



where \l) is the intermediate electronic state with en- 
ergy £; and S is a damping constant. The current oper- 
ator j ia (k) = ^Y,a( A ai^) p La s ^ + h - c -) describes the 
dipole transition between Mn Is and 4p orbitals. Here, 
Piaa and Si a are the annihilation operators of electron in 
Mn 4p and Mn Is orbitals, respectively, with spin a and 
Cartesian coordinate a. A a (k) is the coupling constant 
between the current and the photon. The contribution 
from the quadrupole transition is small, because the in- 
version symmetry is preserved in this systemEj As shown 
in the next section, the anisotropy of the tensor element 
in the anomalous scattering factor determines the polar- 
ization dependence of the scattering intensity. 

The anomalous part of the atomic scattering fac- 
tor Afi a j3(k', k") is calculated in the microscopic elec- 
tronic model. We consider a MnC>6 octahedron and 
introduce the following orbitals in a Mn ion: {Is, 
3d 7 (7 = 70+,7o_), 4p 7 (7 = x,y,z) }, where |3d 7e+ ) = 
cos(6K t 72)|3cZ 32 2_ r 2)+sm(6>( t V2)|3d 2; 2_ a 2) and |3d 78 _) = 
-sin(6»(*V2)|3d 3z 2_ r 2) + cos{0^ /2)\3d x 2_ y 2) . Six O 2p 
orbitals, which contribute to the c-bond with the Mn 
3d orbitals, are also considered. These are denoted 
as 2pi (i = 1 ~ 6) where 2p 1(2 ), 2p 3(4 ) and 2p 5(6 ) 
are 2p x orbital at the position of (+(— )a/2, 0, 0), 2p y 
at (0, +(-)a/2,0), and 2p z at (0, 0, +(-)a/2), respec- 
tively, with Mn-0 bond length a/2. The position of the 
Mn ion is chosen to be (0,0,0). By utilizing the irre- 
ducible representation in Oh group, these O 2p orbitals 
are recombined as {2p 7e+ , 2p 7e _ , 2p x , 2p y , 2p z , 2p r 2}. For 
example, 2p 3z 2_ r 2 and 2p x orbitals are represented as 
2p 3z 2_ r 2 = -^(i (_2 Pl + 2p 2 - 2p 3 + 2p 4 ) + 2p 5 - 2p 6 ) , 

and 2p x = -j^(2pi + 2p 2 ) ■ 

On these bases, we set up the following Hamiltonian: 



H — H n + H t + H 3 d-4 P + H4 P -2 P + H c 



3d-3d 



(15) 



The first and second terms describe the energy level in 
each orbital and the electron transfer between Mn 3d and 
O 2p orbitals, respectively, and are given by 



H + H t = ^2 e d d\ je<J d le<J + ^ e P P\ a P ia 

7fl CT 7 cr 

+ V3t 3d ^ 2 p ^ ( d \e<?P~te<7 + h.C. ] , (16) 



where, d 7e(T , and pro- are the electron annihilation opera- 
tors in Mn 3d, and O 2p orbitals, respectively, with spin 
a and orbital 7e(= 70+, 70-) and 7(= x, y, z). The third 
and fourth terms in Eq.(|l5|), which describe the Coulomb 
interaction between Mn 3d and Mn 4p orbitals and Mn 
4p and O 2p orbitals, respectively, are given by 

H^d-ip + H^p-ip = ^ v (3d le , 4p 7 ) n(3d 7e ) n(4p 7 ) 

78 7 

+ ^V(4p 7 ,2p 7e ) n(4p 7 ) n h {2p le ) . (17) 

78 7 

n(3d 7e ) and n(4p 7 ) are the number operators of Mn 3d 
and Mn Ap electrons, respectively, and n^(2p 7 ) is the 
number operator of the O 2p holes . The explicit forms 
of the Coulomb interaction in Eq. ([l7]) are written by 

F(3d 7e±! 4p 7 ) = F (3d,4p) 

± 4F 2 (3d,4p) cos(0 (t) -I 

and 



V(2p 7e± ,4 Pl ) 



ep' 

-e ± cos 

5 



2tt 



2tt 



(18) 



(19) 



m x — +1, m y = —1, and m z — 0. F n (3d, 4p) is 
the Slater-integral between 3d and 4p electrons, and is 
defined by F (3d,Ap) = F^ (3d, Ap) and F 2 (3d,4p) = 
^.F( 2 )(3d,4p). e = Ze 2 /a and p = (r 4p )/a where Z = 2, 
and (r^p) is the average radius of Mn 4p orbital. Further- 
more, we introduce the interaction between Mn 3d(4p) 
electron and Mn Is core hole as follows; 



Hcore = ^ v ( ls '3d) n(3d l0IJ ) n(ls a >) 

7efcr' 

+ Yl V ( ls ^P) n ^Pl°) n i ls <y') , 



(20) 



where both V(ls,3d) and V(ls,Ap) do not depend on 
the orbitals. The last term in Eq.(in|), H 3 d-3d, is the 
interaction between Mn 3d electronalj defined by 



H 3l i- 



3d-3d 



U Y «(3d 7eT ) n(3d lgl ) 



78 



+ U'n(3d la+ ) n(3d 7e _) 



' d 7fl 1 a 1 d 7 



7e+tr u ' 7e _ C r' u '7e+o- "Ts-er 



2g °78 1 



(21) 



(22) 
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where U, U' and K are the Coulomb and exchange in- 
teractions between 3c? electrons with K = (U — U')/2. 
Jh describes the Hund coupling between the localized 
t2 g spin (St 2 ) with S = 3/2 and e g spin defined by 
s~<e = \ H<j'a" d \ „>{v)o'<y"d le< ," ■ The electron transfer 
{tip-2 P ) between Mn Ap and O 2p is not included in the 
model. As mentioned above, the O 2p orbitals are recom- 
bined by utilizing the irreducible representation in Oh 
group, and the Mn 3c? and Mn Ap orbitals are decoupled 
when t3d-2p and t± v -2p are taken into account. There- 
fore, the essential conclusion about the anisotropy of the 
atomic scattering factor being based on the Coulomb in- 
teractions is not changed, although a large overlap be- 
tween Mn Ap and O 2p orbitals may reduce the Coulomb 
interaction between Mn 3c? and Mn Ap and increase the 
one between Mn Ap and O 2p. 

It is worth noting that H^d-Ap and H^ p -2p provide 
the anisotropy of the atomic scattering factor in or- 
bital ordered states. The essential points in the mech- 
anism of the anisotropy is the fact that these interac- 
tions depend on the occupied Mn 3c? orbitals. When 
an electron occupies 3d^ z 2_ r 2 orbital (9^> — 0), the 
energy of Ap z orbital is higher than that of Ap x r y \ or- 
bital by 6F 2 (3c?, Ap) due to V(3<£y e ,4p 7 ), as a result, 
the scattering intensity near the K-edge is dominated by 
4px(y) orbital. Furthermore, \3d^ z2 _ r 2) state is strongly 
mixed with \3d\ z2 _ r2 3dl, 2 _ y2 p x 2_ y2 ), where \p x 2_ y 2) de- 
scribes the state where a hole occupies the O 2p x 2_ y2 
orbital. Through the inter-atomic Coulomb interaction 
V(2p le) Ap 1 ), the energy of Ap z orbital becomes higher 
than that of Ap x ^ orbital by j^sp- Therefore, the two 
interactions due to the intra- and inter-site Coulomb in- 
teractions result cooperatively in the anisotropy of the 
scattering factor. In the following numerical calculation, 
Hip~2 P is often neglected, since the inter-site interac- 
tion is well screened in comparison with the intra-site 
Coulomb interaction. 

The lattice distortion in the MnOg octahedron also be- 
comes the origin of the anisotropy of the scattering factor. 
When there is the Jahn-Teller(JT)-type lattice distortion 
characterized by the difference of the bond lengths in the 
a6-plane and the c-direction, i.e., 5a = a z — a x i y \, the 
energy difference of the Mn Ap x ^ and Ap z is brought 
about. In the case of a z > a x ( y ) 1 the energy of the 
Mn Ap z orbital is relatively stabilized in comparison with 
that of Mn Ap x (y) orbital. Therefore, this contribution 
competes with the above originated from the intra- and 
inter-atomic Coulomb interactions. In LaMnOa where 
the 5a/ a is about 15%j£3 the contributions from the 
Coulomb interactions seems to be weakened by that 
from the lattice distortion. However, the notable lat- 
tice distortion fn-|Mn06 octahedron is not observed in 
Lao.5Sr1.5MnO4.E3 Furthermore, it was recently reported 
that in Lao.ssSro.^MnOa the X-ray scattering intensity 
attributed to the anomalous scattering is observed at 
(030) reflection in the O* phase,c3 where the six Mn-0 
bonds in the MnOe octahedron are equivalent O There- 



fore, we conclude that the Coulomb interactions domi- 
nate the origin of the anisotropy of the anomalous scat- 
tering factor in these compounds and do not introduce 
the contribution from the lattice distortion in the Hamil- 
tonian in Eq.([T5|). 

By adopting the above Hamiltonian, we calculate the 
atomic scattering factor in the configuration interaction 
method. |3c? 7e+ ) and |3c? 7()+ 3d 7() _2p 7() _) states are intro- 
duced as the bases of the initial and final states, and 
|3c? 7e+ 4p 7 Is) and |3c? 7e+ 3c? 7() _2p 7£) _4p 7 Is) states are as 
the bases of the intermediate states. The second term 
in the right hand side in Eq. ( O ) is only considered. 
The parameter values in the Hamiltonian are chosen to 
be e P - e d = 13.0ey, e p - e d = 2.5eV, t 3 d-2 P = 1.2eV, 
V(3d, Is) = V(Ap, Is) = 7.0eV : U = 8.5 e y, U' = 6.5eV 
and Jh = l.OeV. The Slater-integrals between Mn Ap 
and Mn 3c? electrons are chosen to be F (3c?, Ap) = 5.0eV 
and F 2 (3o!, Ap) = 0.3eV^ which are the same order as 
ones between Mn 3c? electrons JF n (3d, 3d) evaluated by 
the photoemission experiments.Ej The damping parame- 
ter 5 in Eq. (jl^) is taken to be 0.5eV. 

In Fig. 2, we present the energy dependence of the real 
and imaginary parts of the anomalous scattering factor. 
The edge of the lowest peak of Af" aa corresponds to the 
Mn K-edge. Near the edge in Fig. 2(a), where 3d 3z 2_ r 2 
orbital is occupied, the anisotropy is clearly shown. The 
energy position of &-f" x ,s is lower and its weight is 
larger than those of Af" zz . In Fig. 2(b) where 3d x 2_ y2 
orbital is occupied, the anisotropy in the main peak is 
entirely opposite to that in Fig. 2(a). In both Figs. 2 
(a) and (b), the spectra of Afi XX and &fi yy are identical 
as expected. 



IV. POLARIZATION DEPENDENCE OF THE 
SCATTERING INTENSITY 

By using the general formulae of the scattering 
intensity (Eq.(Q)-(ra)) and the atomic scattering fac- 
tor calculated in the microscopic electronic structure, 
we numerically calculate the polarization dependence 
of the normalized scattering intensity I(tp,tpA) — 
I(ip,ipA)/(N] iL e 4: /(mc 2 ) 2 \FA\ 2 ), where N A is the number 
of the Mn atom. The energy of the incident photon beam 
is fixed at that of the K-edge in a Mn 3+ ion. In our model, 
it is chosen to be the energy at which the lowest compo- 
nent of Af" aa (a — x,y,z) has a half value of its maxi- 
mum. The normal part of the scattering factor of a Mn 3+ 
ion is chosen as foim/\A aa \ 2 = 5eV~ x which is about 5 
times larger pthan the typical value of |A/j aa |m/|j4 aa | 2 
at the edge.Eil The following part in this section is di- 
vided into three subsections and the numerical results 
calculated in several orbital ordered cases are presented. 
The qualitative features in the following numerical re- 
sults, e.g. periodicity and phase of the oscillations as a 
function of tp, are independent of details of the micro- 
scopic calculation shown in the previous section and are 
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related to symmetry of the orbital ordering and the ex- 
perimental arrangement, although the detail quantitative 
ones depend on the adopted interactions and parameter 
values in the microscopic model. 



A. ferro-type orbital ordering in the simple cubic 
lattice 



We consider the ferro-type orbital ordering (orbital-F) 
in the simple cubic lattice where a kind of orbital is oc- 
cupied in each Mn site and demonstrate how to identify 
types of orbital through the polarization analyses. The 
numerical results shown in this subsection will help us 
to understand the polarization dependence of the funda- 
mental reflections in the antiferro-type orbital ordering 
shown in the next subsection. 

In orbital-F case, the structure factor per unit cell 
(Eq. (0)) is given by jfF al3 = f Aa /3 at the (hkl) = 
(n x n y riz) reflection with integer n,i(i = x,y,z). Here, 
(hkl) is the indices in the cubic coordinate defined as 
K = h—a + k—b + l—c where a is the lattice pa- 

a a a 1 

rameter in the cubic cell. fAa/3 is the sum of the nor- 
mal and anomalous parts in a Mn 3+ ion with orbital A; 
SapfoA + AfAa/3- In Fig. 3(a), we present the scatter- 
ing intensity with 9^ = which means that an electron 
occupies the 3d 3z 2_ r 2 orbital in each Mn 3+ ion. It is cal- 
culated in the cases of (ei, 62,63) = (b, c, a) (azimuthal 
scan around a-axis). In the case of (ei, e2, e^) = (c, a, b) 
(azimuthal scan around 6-axis), ip dependence is entirely 
opposite to that in Fig. 3(a). On the contrary, the inten- 
sity is independent of p in the case where the azimuthal 
scan is performed around c-axis. We first consider the 
contribution from the normal part of the scattering fac- 
tor. The polarization dependence of the intensity is rep- 
resented by I (ip, (pa) = | foA cos <Pa | 2 , which is indepen- 
dent of <p due to the scalar character of f A- For con- 
venience, the value of ipA = is defined Iq = |/oa| 2 - 
This pa dependence is interpreted that the scattering at 
Pa — (tt/2) is accompanied with the polarization of 
(A', A") = (a,a)((a,nj), where A' (A") describes the po- 
larization of the incident (scattered ) photon. Therefore, 
the polarization dependence in Fig. 3(a) is dominated 
by the normal part of the scattering factor and slightly 
modified by the anomalous part. 

Let us examine the scattering intensity at pa = in 
Fig. 3(a). This is represented by 



/(^,^ = 0) = |/oa + A/| 5 



with 



A/ = Af Ayy cos 2 <p + Af A 



, sin 2 (p 



(23) 



(24) 



Since the normal part of the scattering factor dominates 
the scattering intensity, it is approximated as I(p, pa — 
0) ~ foA + where Af' A is negative at the 



edge. As a result, the real part of the anomalous com- 
ponent decreases the scattering intensity from Iq. In the 
case where 3d 3z 2_ r 2 orbital is occupied, |A/^. 



Axx(yy) 



has 



a larger value at the edge in comparison with |A/^ ZZ | 
(see Figs. 2(a) and (c)). Therefore, the reduction of the 
scattering intensity from Iq is remarkable around ip = 
and 7r, where the electric vector in the incident photon is 
parallel to the 6-axis. 

In Fig. 3(b), we present the polarization dependence of 
the scattering intensity in several orbital-F cases at pa = 
0. With increasing 9^ from 9 A — 0, the ip dependence 
becomes weak and disappears at 9 A = 2w/6 (3d y 2_ z 2). 

With further increasing 9 A , the phase of a modulation 
becomes inverse, since Af Azz dominates the scattering 
factor near the edge. The inversion of a modulation in 
I(p,PA — 0) occurs at 9 A — An/6 and in the cases 
of (ei,e2,e*3) = (c,a,b) and (a, b, c), respectively. By 
combining the azimuthal scan around a-, b- and c-axes, it 
is possible to identify the occupied orbital in the orbital-F 



B. antiferro-type orbital ordering in the simple 
cubic lattice 

Next, we consider the antiferro-type orbital ordering 
(orbital- AF), where two kinds of orbital sublattice, de- 
noted by A and B, exist. In the simple cubic lattice, 
there are three kinds of orbital-AF, that is, layer-type 
(orbital- A AF), rod- type (orbital-CAF) and NaCl-type 
(orbital-GAF). These notations are defined— by, analogy 
with types of antiferromagnetic structures .EZrCJ In each 
case, two kinds of the reflection point exist, termed the 
fundamental reflection at (hkl) = (n x n y n z ) and the or- 
bital superlattice reflection at (hkl) = (n x n y n z + h) 
for orbital-CAF, (n x + | ny + j n z) for orbital- AAF and 
(n x + 2 n y + 5 n z + |) for orbital-GAF. By examining 
the superlattice reflection, the above three types of the 
orbital-AF is able to be distinguished. The structure fac- 
tor per unit cell is written as -kF a p — (fAa/3 + /bq/s) for 
the fundamental reflection and jjF a p = (fAa/3 — fsa/s) 
for the superlattice reflection. For convenience, we define 

f+a/3 = \(fAaf3 + f Baf3) and f- a /3 = \(f Aaj3 ~ fBap): 

and term the former and latter the ferro- and antiferro- 
components of the scattering factor, respectively. We ne- 
glect the difference of the normal part of the scattering 
factors of Mn 3+ ions between the A and B sublattices. 
Therefore, at the orbital superlattice reflection, only the 
antiferro-component of the anomalous part contributes. 

In Fig. 4(a), the polarization dependence of the 
fundamental reflection in the orbital-AF case is pre- 
sented. The orbital state in the two sublattices is 
{®A ->@b) ~ (47r/6, — 47r/6) which corresponds to the 
(3d 3 j,2_ r 2 , 3d 3;r 2_ I ,2)-type orbital ordering. The coordi- 
nate in the laboratory system is chosen to be (ex , 62, 63) = 
(6, c, a). The explicit formula of the intensity in Fig. 4(a) 



5 



at (fiA — is given by Eq. (^3|) with 

Af = Af +yy cos 2 Lp + Af +zz sin 2 <p . (25) 

We note that A/a vv and A/a zz in Eq.(p^) are replaced 
by the ferro-components of the scattering factor, i.e., 
Af +yy and Af +zz in Eq.(p5|), respectively. The reduction 
of the intensity from Iq becomes remarkable around ip = 
tt/2 and 3ir/2 where the electric vector in the incident 
beam is parallel to the c-axis. This is because the ferro- 
component of the scattering factor ^(AfAaa + AfBaa) 
is larger for a — z in comparison with that for a = x{y) 1 
since the 3(^2 _ r 2 and 3d 3x 2_ r 2 orbitals are almost elon- 
gated in the afo-plane. When the azimuthal scan is per- 
formed around c-axis, the intensity is almost independent 
of Lp, because the ferro-component of the scattering factor 
along the a- and 6-axes is identical under the condition 

of9® = -0$. 

In Fig. 4(b), the orbital dependence of the intensity at 
the fundamental reflection with ipA = is shown. With 
increasing 9^ from = 7f/6, the ip dependence be- 
comes weak and is almost smeared out at (9^,9^) = 
(7r/2, — 7r/2). With further increasing 9j±, the phase of 

the modulation is reversed. In the case of 9^ > tt/2, the 
ferro-component of the scattering factor in the c-direction 
becomes dominant in comparison with that in the ab- 
plane. We conclude that through the measurement of 
the phase of the polarization dependence, we can deter- 
mine whether 9^ is larger or smaller than tt/2. 

In Fig. 5(a), the orbital superlattice reflection, which 
appears only in orbital- AF ordering is shown. The orbital 

state is (9^ , 9$) — (4tt/6, — Att/6) which corresponds to 
(3d 3y 2„ r 2, 3d 3x 2_ r 2) ordering. The coordinate in the lab- 
oratory system is chosen to be (ei, €2, eb) = (c, ^75 (a — 
b), -^=(a + b)). The azimuthal scan in this coordinate was 
carried out in Lao.5Sro.5Mn04 and LaMn03. 

EE 

The ex- 
plicit formulae of the polarization dependence at lpa = 
and 7r/2 are given by 

I(tp, Lp A = 0) = I Af_ zz cos 2 tp 

+ i(A/_ M + A/_ ro )sin 2 vf , (26) 

and 

I(ip,Lp A = tt/2) = I ^ ( Af— X x + Af- yy ) sin^costf 

+ -{Af_xx + A/_ TO - 2Af_ zz ) 

x sin tp cos Lp sin 9 | 2 , (27) 

respectively. The intensity is expressed only by the 
antiferro-component of the scattering factor. At this re- 
flection, the information of the orbital ordering is derived 
without disturbance of the normal part of the scatter- 
ing factor, in contrast with the case of the fundamen- 
tal reflection. When the orbital ordering is assumed to 



be [0$,eg = -6$), the condition, Af Axx = Af Byy , 
MAyy = AJbxx and Af Az z = Af Bzz , is satisfied. In 
this case, I(Lp,Lp A = 0) in Eq.(p6|) becomes zero, and 
IJ&VA =V 2 ) = \\(Af Ax x- Af A yy) sin Lpcos9\ 2 in Eq. 
(P7|). These expression gives a square of the sinusoidal 
symmetry in the polarization dependence as shown in 
Fig. 5(a). By utilizing this characteristic feature of the 
polarization dependence, we can judge whether the con- 
dition {9 9^ = is satisfied or not. For example, 
as shown in Fig. 6, when we assume different types of 
the orbital- AF ordering {9^, 9$ = 9^ + tt), the polar- 
ization dependence shows a quite different feature from 
that in Fig. 5. 

In Fig. 5(b), the polarization dependence at LpA = tt/2 
in several orbital-AF cases is presented. Here, the con- 
dition {9^a i^b = — @a) i s assumed. The scatter- 
ing intensity changes with 9y and becomes maximum 

around 0y = tt/2. The polarization dependence is, 
however, given by a square of the sinusoidal function 
in all cases. Ip^the azimuthal analyses performed in 
Lao.sSri.sMnOj 2 -!!, the experimental data are well fitted 
by the above function. Therefore, we conclude that the 
orbital ordering in Lao.sSri.sMnO^ satisfies the condi- 
tion (9a, 9^ = —Ojl), i-e., the A and B orbitals are 
symmetric with respect to replacement of x(y) by y(x). 
However, we can not determine which type of the orbital 
ordering, (3d 3x 2 _ r 2 , 3d^ y 2 _ r 2) or (3d z 2_ x 2 i 3d y 2_ z 2) is re- 
alized, through the measurement of this reflection. 



C. orbital ordering in charge ordered state 

The charge ordering, where Mn 3+ and Mn 4+ 

ions are alternately aligned in the Mn02 plane, have 
been observed ia|JjLi-|- x A x MnC)3 near x = 0.5 and 
Lao. 5 Sr 1 . 5 Mn04.E3ii3'E3 In the latter .compound, it is ac- 
companied with the orbital ordering^ and the CE-type 
spin ordering. In this subsection, we investigate the po- 
larization dependence of the scattering intensity in the 
charge ordered phase and propose the possibility to iden- 
tify the orbital ordering in this phase. 

We consider the alternating alignment of Mn 3+ and 
Mn 4+ ions in the MnC>2 plane and introduce two kinds 
of the orbital sublattice for Mn 3+ ion, denoted by A and 
B. It is assumed that a Mn 3+ ion has A(B) orbital, when 
the spin in a Mn 3+ ion is parallel to lhat.in.jta neighbor- 
ing Mn 4+ ions in the a(6)-direction.EOl3'BEl This type 
of orbital ordering is consistent with the CE-type spin 
structure. Here, we pay our attention to the reflection at 
(h k I) — (n x + 5 n y + \ n z) termed the charge-order re- 
flection, which appears due to the existence of the charge 
ordering. The structure factor at the reflection is given 
by 

F a j3 = N(fAaf3 + fBafi ~ 2fi a p) , (28) 



G 



where fA(B)a/3 1S the scattering factor defined in the or- 
thohombic coordinate as follows 

(2~{flxx ~\~ flyy) ^(.flxx flyy) ^ \ 
\{flxx-flyy) \{flxx+flyy) , (29) 

" f lzz J 

for I — A and B. is the atomic scattering factor of 
Mn 4+ ions. We emphasize that the difference of the nor- 
mal part of the scattering factor f A + foB~ 2/o4 is of the 
order of (one electron) /(Mn atom), which is much smaller 
than the typical value of the anomalous part at the edge. 
Furthermore, the anomalous part of the scattering factor 
of a Mn 4+ ion is scalar. Its value is smaller than that 
of Mn 3+ ions near the K-edge of a Mn 3+ jon, because 
the edge is lower than that of a Mn 4+ ionO Therefore, 
by utilizing this reflection, we can obtain the information 
of the ferro-component of the scattering factor without 
disturbance of the normal part. 

In Fig. 7(a), the polarization dependence of the 
scattering intensity at the charge-order reflection is 
presented. The types of the orbital is (6^,6g) — 
(47r/6, — 47r/6). The coordinate in the laboratory system 
is chosen to be (ei, £2, £3) = (c, -j=(a — &), -j=(a + b)). We 
neglect the contribution from Jqa + Job — 2/o4, and as- 
sume that the anomalous part of the scattering factor in 
a Mn 4+ ion as a scalar, that is, Afj± a p = Af^dap, and the 
value of A/ 4 is a half of the maximum value of |A/;Q, a | 
(/ = A, B) (a = x,y,z) in a Mn 3+ ion.El The explicit 
formula at ip A — is given by 

I{ip,<fA = 0) = — \(Af +xx + Af +yy )sm 2 ip 

+ 2A/ +M cos 2 <^-2A/ 4 | 2 . (30) 

In the case of <p — and 7r in Fig. 7(a), the electric 
vector in the incident beam is parallel to the c-axis. It 
is worth noting that the polarization dependence in Fig. 
7(a) is more remarkable in contrast with that shown in 
Fig. 4(a), because the scalar component of the scattering 
factor is smaller in the present case. Furthermore, the in- 
tensity in the case where the electric vector is parallel to 
the afo-plane becomes maximum, that is, the polarization 
dependence is opposite to the case in Fig. 4(a). This is 
attributed to the fact that the polarization dependence is 
dominated by the anomalous term of Mn 3+ ions. On the 
other hand, in Fig. 4(a) the interference term between 
the normal and anomalous parts determines the polar- 
ization dependence. As shown in Fig. 7(b), the phase of 

the polarization dependence for 0$ > tt/2 is opposite to 

that for 9jp < n/2. It is concluded that in the charge 
ordered state, by utilizing the polarization dependence 
of the charge-order reflection the ferro-component of the 
scattering factor in orbital-AF state is identified easier 
than that in the case of fundamental reflection. 



V. SUMMARY AND DISCUSSION 

In this paper, we have theoretically investigated the 
anomalous X-ray scattering as a probe to detect the 
orbital ordering in manganites. Through the polar- 
ization dependence of the scattering factor, we showed 
how to identify several types of the orbital ordering. 
In particular, we paid our attention to the method 
to distinguish the two types of the orbital orderings; 
(3d 3a: 2_ r 2, 3d 3y 2_ r 2) or (3d x a_ x a, M y 2_ z 2). 

At the superlattice reflection due to the antiferro-type 
orbital ordering with the condition 0g = — 0^ \ the 
azimuthal angle dependence is represented by a square 
of the sinusoidal function and its intensity depends on 
types of orbital (Fig. 5). We found that it is difficult 
to determine which orbital ordering (3d 3x 2_ r 2,3d 3y 2_ r 2) 
or (3d z 2_ a .2, 3d y 2_ z 2) is realized by analyzing the results 
at the reflection. This difficulty is attributed to the fact 
that the antiferro-component h(AfA a ~ Afs a 0) of the 
scattering factor is observed at the superlattice reflection. 
Then, we have proposed two kinds of polarization anal- 
yses where the ferro-component of the scattering factor, 
i.e., \{Af Aa i3 + A f Baft), is derived. 

The first one is the polarization analyses at the fun- 
damental reflection (Fig. 4), where the scattering factor 
is represented by sum of the normal part and the ferro- 
component of the anomalous part. The interference be- 
tween them gives rise to the polarization dependence of 
the scattering intensity. Since the normal part is much 
larger than the anomalous one, the orbital ordering re- 
flects on the modulation in the polarization dependence. 
In order to detect the modulation, the reflection plane is 
required not to include the heavy ions, such as a La ion, 
because the large values of the normal part disturbs to 
detect the modulation. The reflection with the large scat- 
tering angle characterized by the large number of (hkl) 
is also suitable, because the normal part of the scattering 
factor is reduced with increasing the scattering momen- 
tum, in contrast with the anomalous part. 

The second proposal to obtain the ferro-component of 
the scattering factor is the analyses at charge-order re- 
flection points in charge ordered phases (Fig. 7). At the 
reflection, the normal part is almost canceled out. The 
scattering factor is dominated by the ferro-component of 
the anomalous scattering factor in a Mn 3+ ion and the 
anomalous part in a Mn 4+ ion . The latter is a scalar 
and its value is smaller than that of a Mn 3+ ion near 
the Mn 3+ K-edge. Therefore, by analyzing the polariza- 
tion dependence of this reflection, the ferro-component 
of the scattering factor is observed. It is noted that the 
basic considerations in the above theoretical proposals 
are independent of details in the numerical microscopic 
calculation and are based on symmetry of types of the 
orbital ordering and the experimental arrangement. 

In conclusion, we have shown that in the anomalous 
X-ray scattering technique, types of the orbital order is 
able to be identified by selecting the adequate reflection 
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point and analyzing the polarization dependence of the 
scattering intensity. 
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Figure captions. 

Fig. 1: 'ftp - ! schematic view of the experimental 
arrangementEHtS which consists of a sample crystal (S), 
analyzer crystal (A) and photon detector (D). The polar- 
ization dependence of the scattering intensity is measured 
as a function of the azimuthal angle (<p) and the analyzer 
angle (</?a). k' (k") is the photon momentum and it' (n") 
and a' (cr") are the polarization of the incident (scattered) 
photon. 

Fig. 2: The real and imaginary parts of the scattering fac- 
tor ((Afi) xx ( yytZZ )m/\A\ 2 ) in the cases where the following 
orbitals are occupied: (a), (b) = (3d 3z 2_ r 2), and (c), 
(d) 9^' = 7r (3d x 2_ y 2). The straight and broken lines show 
(Afi) xx ( yy \ and (A/i) zz , respectively. The parameter value 
of 8 in Eq.(p^) is chosen to be 0.5eV. The origin of the en- 
ergy is taken to be arbitrary. The spectra of (Afi) xx and 
(Afi) yy are identical. 

Fig. 3: The polarization dependence of the scattering inten- 
sity in the orbital-F case. The coordinate in the laboratory 
system is chosen to be (ei, e 2 , ez) — (b, c, a), where (a, b, c) 
represents the cubic crystallographic coordinate, (a): The 
orbital is chosen to be 9^ — (3d 3z 2_ r 2). (b) : The az- 
imuthal angle dependence at ifA — in several orbital-F 
cases. 

Fig. 4: The polarization dependence of the scattering in- 
tensity of the fundamental reflection in the orbital-AF 
case. The coordinate in the laboratory system is chosen 
to be (ei, e 2 , £3) = (b, c, a) where (a, b, c) represents the cu- 
bic crystallographic coordinate, (a) : (6^ = 4-7r/6,6l^ = 
— 47r/6) orbital ((3d 3y 2_ r 2, 3d 3x 2_ r 2)) is assumed, (b): The 
azimuthal angle dependence at ipA = in several orbital- 
AF cases with (ffl ,0 { g = -Of). 

Fig. 5: The polarization dependence of the scattering inten- 
sity of the orbital superlattice reflection in the orbital-AF 
case. The coordinate in the laboratory system is chosen to 
be (ei, e 2 , e 3 ) = (c, -^{a—b), -±=(a + b)), where (a, b, c) rep- 
resents the cubic crystallographic coordinate, (a): (9^ — 
47r/6, O^g = — 47r/6) orbital ((3d 3y 2_ r 2, 3d 3x 2_ r 2)) is as- 
sumed, (b) : The azimuthal angle dependence at ipA = tt/2 
in several orbital-AF cases with (0^,8^ — —9 a)- 

Fig. 6: The polarization dependence of the scattering in- 
tensity of the orbital superlattice reflection in the orbital- 
AF case. The coordinate in the laboratory system is cho- 
sen to be (ei,e 2 ,e 3 ) = (c, ^=(a — 6), ^=(a + &)), where 

(a, b, c) represents the cubic crystallographic coordinate, 
(a): {9a = Ot"B — 7r ) orbital ((3d 3z 2_ r 2, 3d x 2_ y 2)) is as- 
sumed, (b) : The azimuthal angle dependence at ip a = tt/2 
in several orbital-AF cases with {9f,9$ = 9f +7r). 

Fig. 7: The polarization dependence of the scattering inten- 
sity at the (n^-l-i n y + ^ n z ) reflection (charge-order reflec- 
tion) . The coordinate in the laboratory system is chosen to 
be (ci,e 2 , e 3 ) = (c, -^j(a-b), -±=(a + b)) where (a, b, c) rep- 
resents the cubic crystallographic coordinate, (a): (9^ — 
4n/6,9^ = — 47r/6) orbital ((3d 3y 2_ r 2, 3d 3x 2_ r 2)) is as- 



sumed, (b) : The azimuthal angle dependence at if a = 
in several orbital-AF cases with (9^,9^ = —8^). 
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Fig. 1 S.Ishihara et al. 
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Fig. 2 S.Ishihara et al. 
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Fig. 3 S.Ishihara et al. 
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Fig. 4 S.Ishihara et al. 
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Fig. 5 S.Ishihara et al. 
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Fig. 6 S.Ishihara et al. 



(a) e (t) A =W6 (3dj v ' 3d 3xV ) 




0.000 j- 1 1 1 1 1 1 — 1 — 1 

71/2 71 37C/2 In 



9 



Fig. 7 S.Ishihara et al. 



